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Abstract. 

We continue the study of the entanglement entropy of two disjoint intervals 
in conformal field theories that we started in J. Stat. Mech. (2009) P11001. We 
compute Trp^ for any integer n for the Ising universality class and the final 
result is expressed as a sum of Riemann-Siegel theta functions. These predictions 
are checked against existing numerical data. We provide a systematic method 
that gives the full asymptotic expansion of the scaling function for small four- 
point ratio (i.e. short intervals). These formulas are compared with the direct 
expansion of the full results for free compactified boson and Ising model. We 
finally provide the analytic continuation of the first term in this expansion in a 
completely analytic form. 
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1. Introduction 

In this manuscript we continue the study of the entanglement entropies of two disjoint 
intervals in conformal field theories (CFT) that we started in Ref. pQ. Let us consider 
a one-dimensional critical system whose scaling limit is described by a CFT of central 
charge c, and a partition into two parts A and its complement B. We define the 
reduced density matrix pa — Tr^p, where p is the density matrix of the entire system. 
The entanglement between A and B is characterized by the moments of pa, i-e. Tr//|, 
from which through the replica trick [2J|3] one obtains the Von Neumann entanglement 
entropy [3] 



d 

Sa = -TrpA log pa = - -q^Ttpa 



(1) 



For integer n, Trp^ is proportional to the partition function on an n-sheeted Riemann 
surface with branch cuts along the subsystem A, i.e. Trp^ = Z n {A)/Zi where Z n {A) 
is the partition function of the CFT on a conifold where n copies of the manifold 
1Z = system x R 1 (or system x S 1 at finite temperature) are coupled along branch 
cuts along each connected piece of A at a time-slice t = 0. 

In the case when A is a single interval of length I in an infinite system, the n- 
sheeted Riemann surface is topologically equivalent to the complex plane on which 
it can be mapped by a uniformizing mapping. By studying the transformation 
of the stress energy tensor (i.e. the response of the partition function to a scale 
transformation) under this uniformizing map, one obtains [3J |3] 

Trp% = c n r^ n - 1 ' n \ (2) 

from which the replica trick ([I]) leads to the famous result (c/3) log^ [5][S] and to the 
the full spectrum of pa [3- 

When the subsystem A consists of N disjoint intervals (always in an infinite 
system), the entanglement entropy is sensitive to the full operator content of the 
theory and not only to the central charge c. The n-sheeted Riemann surface 7Z n ,N 
has genus (n — 1)(N — 1) and cannot be mapped to the complex plane so that the 
CFT calculations become much more complicated. Only few results are analytically 
available. For two intervals, i.e. A = A\ U A2 = [ui,vi] U [^2,^2] (without loss of 
generality we assume V2 > U2 > v\ > u\), global conformal invariance leads to the 
scaling form 

^ - % = 4 ( 7 JT'^'^I ^ § ( "" Al) > (3) 



(4) 



Z i "V ( v i - ui)(v 2 - u 2 )(v 2 - ui)(u 2 - V\) 

where x is the four-point ratio (for real Uj and Vj, x is real) 

(v-l - ui)(u 2 - u 2 ) 
x = 

(U 2 - Ul)(v 2 - Vi) ' 

In Ref. pQ we showed that for a boson compactified on a circle of radius R (a Luttinger 
liquid field theory), the universal scaling function J- n {x) is 

_ e(o| n r)e(o|r/ n ) 
n[) ~ [e(o|r)p ' [) 

where T is an (n — 1) x (n — 1) matrix with elements 

(6) 



2i / k . 

Trs = — 2^ Sm [ ) Pk/n COS 



k = l 



II 



2tt — (r — s) 
n 
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and 



•v - 



gjjXg/, 1 - y; !; 1 - x ) 

2 F 1 (y, 1-y; l;x) 



(7) 



r\ is a universal critical exponent proportional to the square of the compactification 
radius R 1 while is the Riemann-Siegel theta function 



with z a generic complex vector with n— 1 components (we denote vectors in bold letter 
to easily distinguish them from scalars). For n = 2, .7-2 (x) reduces to the result in [8]. 
This is the only result for F n (x) known in full generality, and many other important 
partial results are known [l^IiniinjinillllllEa^ both 
analytically and numerically. 

However, although the moments of the reduced density matrix have been obtained 
for all integer n, the analytic continuation to complex n is still beyond our knowledge 
and so is the Von Neumann entropy. Some results have been obtained in Ref. |I] from 
Eq. ([5]) in the limit of small r\ and of small x (the latter only numerically, except for 
few specific values of rj), and the two limits do not commute. 

In this manuscript we extend the known results in several directions. We provide 
a systematic method that gives the full asymptotic expansion of T n (x) for small x 
(i.e. short intervals). This method is valid for any conformal invariant theory (even 
in the presence of boundaries and at finite temperature) and the obtained results are 
tested against known results (Sec. [3]). Furthermore, it provides the coefficients of any 
expansion when the lengths of the intervals are smaller than any other length of the 
system, such as finite size, inverse temperature etc. In these cases the entanglement 
entropy is sensitive to all the CFT data on the appropriate Riemann surface. We 
derive the scaling function T n {x) for all n for the Ising universality class, generalizing 
the result in Ref. [T5] for n = 2 (in Sec. [4]). We extend the small x expansion to the 
second non vanishing order both for the compactified boson and for the Ising model 
(Sec. [5]). We derive a full analytic formula for the first non vanishing term in the small 
x expansion of the Von Neumann entropy (Sec.|6|. 

2. Summary of results 

This section contains a summary of our main results, that are derived in the following. 
2.1. Short-length expansion 

We develop a general framework to calculate the short-length expansion of Tr//| that is 
very similar to the standard short-distance expansion for correlators. This formalism 
(that extends and generalizes a recent approach by Headrick [TH]) applies every time 
in which A consists of one or more intervals {I p } whose lengths £ p are much smaller 
than all the other characteristic lengths ~ L (for example, the separation between 
the intervals, the size of the system or the inverse temperature). The short length 
expansion provides a series of powers of £ p /L whose terms are universal and encode 
all the data of the CFT such as the scaling dimensions and the operator product 
expansion (OPE) coefficients. This expansion makes manifest the common belief that 
the entanglement entropies of disjoint intervals encode information about the full data 




exp [ 




(8) 



mEZ- 1 



of the CFT. 
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When these arguments are specialized to the case of two intervals (of lengths i\ 
and £2 whose centers are at a distance r) in an infinite system, we obtain 

Ti 'Pa = ,„ wfif„._i/„) E '«,>(»)• 0) 

{kj} V 7 



(44) c/6(r 



where the sum is over all the set of possible operators of the CFT with 

conformal dimensions (A^. , A^.). s^.y(n) are calculable coefficients depending on 
the correlation functions of these operators. This is easily converted in a systematic 
small x expansion of the function F n (x). Each set of operators {kj} gives a 

leading contribution at order x^-^ A:,+Aj '. Furthermore, in T n (x) there are analytic 
contributions in x, which do not correspond to any operator and they come expanding 
the prefactor in Eq. ([3| in powers of x. 

In the case of compactified free boson and Ising model, we can expand the exact 
results for any x to obtain 



F n (x) = 1 



( x \ 2a 1 \ 
1 1 54(71) 



(10) 



V4n 2 / 

where a — min[?7, 1/77] for boson and a = 1/4 for Ising. The coefficients S2{n) and 
84(11) are explicitly calculated in the asymptotic expansion and they are identified with 
the contributions of the short-length expansion coming from the two- and four-point 
functions of the most relevant operator in each theory. Contributions from three-point 
function are subleading. s 2 (n) has a simple zero at n — 1, while s^n) has simple zeros 
at n = 1,2,3. In general s m (n) (coming from an m-point function) will have simple 
zeros at any integer n < m. 



2.2. Result for the Ising universality class 

For the Ising model, the scaling function F n (x) is 



T n (x) 



1 



L e(o|r) 



E 

e,S 



e 



(o|r) 



Here 8 is the Ricmann thcta function with characteristic defined as 



e 



£ 

6 



(z|r) 



E 

m£Z' 



exp 



i7r(m + ef ■ T ■ (m + e) + 2?ri (m + ef ■ (z + 5) 



(11) 



(12) 



where z G C" _1 and T is the same as in Eq. e, S are vector with entries and 
1/2. The sum in (e, 6) in Eq. ( fllj ) is intended over all the 2™ _1 vectors e and d 
with these entries. We show that T n (x) = T n (l — x) and that it reproduces known 
analytical and numerical results [T51 [T7] . 



2.3. Analytic continuation of the leading term in the small x expansion 



From the short length expansions (and from the asymptotic expansion of known 
results), the leading term of the small x expansion in T n (x) is 

n-l 

(13) 



and it comes from the two-point function of the most relevant operator of the CFT. 
The integer counts the number of inequivalent correlation functions giving the same 
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contribution (coming most often from degenerate operators). We have Af = 2 for free 
boson and Af = 1 for the Ising model. 

We provide the analytic continuation to complex n for any value of a whose 
derivative at n = 1 is 

s' 2 (l)=Af ^ T ( a + V , (14) 
2W 4r(a+ |) ' y ' 

that gives the asymptotic expansion of the scaling function of the von Neumann 
entropy 



T VN ( X ) ^ ^ n ( X ) = Af (-) ^ + (15) 



<* 0FT(a + l) 



4T(a+§) 



3. Short interval expansion 



In this section we consider Trp^| in a variety of situations in which A consists of one or 
more intervals {I p } whose lengths l p are much smaller than all the other characteristic 
lengths ~ L (for example, the separation between the intervals, the size of the system 
or the inverse temperature) and show that it can be expanded in a series of powers of 
lp/L whose terms are universal and encode all the data of the CFT such as the scaling 
dimensions and the operator product expansion (OPE) coefficients. Our arguments 
are a generalization of those of Headrick 181 to the case of non-zero conformal spin, 
to the next to leading orders in the OPE and to the case of more than two intervals. 
When these arguments are specialized to the case N = 2 in an infinite system, they 
provide a systematic small x expansion of the function !F n (x). 

The basic idea, familiar from other situations in CFT, is that this ratio of partition 
functions can be thought of as a correlator 

^ }1 = (f[l P )n n , (16) 
1 p=i 

in n copies of the CFT on decoupled manifolds 1Z (notice n as power of TZ and not as a 
subscript). More explicitly, in this correlator each I p may be expanded in a complete 
set of local fields at a given point in the interval, for example its midpoint Zpjf] 

n 

{%} J=l 

where {4>k} denotes a complete set of local fields for a single copy of the CFT, and 
z Pj is the point z p on the jth sheet. 

This statement may be understood most easily within radial quantization. In the 
limit when the lengths l p are small compared with the separation between the intervals 
(and other length scales such as the system size, the distance to any boundary or the 
inverse temperature), we may surround the interval I p by a circle on each sheet and 
consider the state \I p ) induced by coupling the CFTs along the branch cut I p . Since 
at this radius the n sheets are distinct, this state lies in the Hilbert space ® r f = {H.j 

% This may be thought of as a modification of the standard OPE in the plane 

4>a{u)4>biv) = Cabc (« " V)- A "- A »+ A c(u - O-^-^+^c^) , 

c 

to the case when when have non-local twist operators in n copies of the CFT. 
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where T-Lj is the space of the jth CFT in radial quantization about z Pj . Each of these 
is spanned by states \4> k ) in one-to-one correspondence with the scaling operators of 
the CFT. That is, we can write 



{k 4 } 



However, by the operator-state correspondence of CFT, this is equivalent to Eq. (17). 

The main point is now that, like the more conventional short-distance expansion 
of two local operators, the coefficients C{ kj y for a given interval are independent of 
the location and lengths of the other intervals in Eq. (16 1. Therefore they may 
be determined by considering the simplest possible case, that of a single interval 
I = (—1/2,1/2) on the infinite line. In this case the partition function itself is trivial, 
but we can also consider the insertion of other operators Oj on each sheet outside the 
circle of radial quantization. Specifically 

^<n°*>*».* = < j u°i)^ = e n<M°i)°i>^ ■ ( i9 ) 

1 3 'J {* ! 3 

Note the appearance of the ratio of partition functions, reflecting the fact that the 
correlators on the left and right sides are evaluated in different geometries. TZj is the 
jth copy of the complex plane. 

In particular, by choosing the Oj to be a complete set of operators <p kj {zj ) as 
Zj — >• oo, we may use orthogonalitjj§] 

(MWuW) = z- 2A "z- 2Xk S kk , , (20) 
where (A k , A^) are the conformal dimensions of to find that 



C {kj } = 



^ z lim z 2 ^ z 2 ^ ^ <H ^ {z 3 ))k^ , (21) 



where the correlation function on the right is evaluated on the n-sheeted surface 
corresponding to a single interval in an infinite system. Note that this correlation 
function vanishes unless the fusion of Ylj 4>kj includes the conformal block of the 
identity. Also their total conformal spin should be an even integer, since I p is invariant 
under a rotation through n about its midpoint. 

We normalize the single interval partition function according to Eq. ([2]), so that 

C{k]} = ^-(c/6)(«-V«)+£,(A* j+ ^) ^ (22) 

where the d{ kj y are dimcnsionless universal numbers. 



In the case where the (j> kj are all primary, the correlator in Eq. (21 1 is simply 
related to the same correlator in the plane by the uniformizing map 

^™ = {jTf') ■ (23) 

which maps each sheet TZj in a wedge of angle 2ir/n. Infinity on the jth sheet is 
mapped to the nth root of unity e 21 ' 3 '", and f'(z) ~ (l/nz 2 ). Thus, for primary 
operators, 

(n^fe))^a=n/ / fe) Afci ^ Afc ^n^(/fe)))c. (24) 



§ This assumes that the operators <p k are all quasiprimary, that is Li</>/s = 0. 
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So, using Eqs. (|21j>, ([22]), and ([24]), we have 

^ , n 

d {kj }=n-^ {Ak i +Ak *\^h{^ ij,n ))c- (25) 

3=1 

For non-primary operators, derivatives and inhomogeneous terms involving the central 



charge may also arise. Eq. (25) is the key result of this section from which all the 
others follow. 

In the limit t — >■ 0, the sum over the other values of {kj} may be organized 
according to how many of them are non-zero (at least in a unitary CFT with non- 
negative scaling dimensions). 

• The leading term as I — > is given by taking all the = 0, that is fa, = 1, 
the identity operator. 

• The next term comes from taking all the kj = except one, say kj = k. However 
in this case fa cannot be primary, since the 1-point function (fa)c would vanish. 
The most interesting case is when fa is a component of the stress tensor, T or T. 



In that case Eq. (21) becomes 



C T = (:: r o = lim zl(T( Zj ))n n , ■ (26) 

(Note the factor of 2/c because the 2-point function of T is normalized to c/2 
rather than unity.) However [5] 

-a (27) 



so that 



and similarly for d^. 

• The next contribution comes from taking two of the kj to be non-zero. Since the 
product must couple to the identity block this requires the two operators to be in 
the same block, and if they are real and primary they must be the same operator. 
Thus, taking kj t — kj 2 — k, we have, in an obvious notation, 

d (Jlh) = n -2(A fc +AO ( ^( e 2™W« )0fe ( e 2«W«)) c 

-2(A fc +A fc ) 

= 2" ■ (29) 

l e 2iriji/n _ gZmjz/n} k j"g— 2mji/n _ e —2mj 2 /n] k 

Note that we must have j\ ^ j2 . 

• Similarly, if three of the kj are non-zero and correspond to primary fields 

{(j) kl ,4>k 2 ,<l>k 3 ) 

= ^ Al - A2 - A3 ^ 1 ^ 2 ^ 3 (^ 1 (e^)^ 2 (e^ a )^ 3 (e^))c (30) 

rl -A 1 -A 3 -A 3 -A 1 -A 2 -A 3ci23 

[e~^ — e^ Za ] Al+A2_A3 [e ^ — e s^]Ai+a 2 -A3 x C y C ij c permutations 
where we shortened Ak p — A p and C123 = ci Cl k 2 k 3 is an OPE coefficient. 
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• The next term involves the 4-point functions at the nth roots of unity, and so on. 
However, the correlation functions with more than 3 points are different for any 



CFT and so Eq. ( 25 ) cannot be simplified in general. 

Note that for a fixed positive integer n, only the iV-point functions with N < n 
contribute. This means that the coefficients in the short length expansion, as analytic 
functions of n, must have zeroes at positive integer values of n < N. 

We now discuss a number of applications of the short length expansion. 

3.1. Single interval in a finite geometry or near a boundary. 

In this case the 1-point functions ((j^k^TZi can be non vanishing. An interesting 
example is the contribution of the stress tensor. From the above we have 

Trp^ = ^-(c/eXn-i/n) (i + . . . + i( n _ l/n)f{{T) + (T)) + •••). (31) 

(The extra factor of n comes from the sum over j.) We recognize the second factor as 
being proportional to the energy density T t t- In fact we can deduce a rather general 
result: the 0(£ 2 ) correction to the Renyi entropy S n = lntr p\ of a single interval 
in a finite system (or at finite temperature) is given by the local energy density: 

SS n = ^(l + l/n)(£)e 2 . (32) 

This can be checked by comparison with the result in Ref. [5] for a single interval in 
a system of length L with periodic boundary conditions 



S n — 







L 


. nf 








sin — 




nj 


7T 


L 



(33) 



by expanding to 0((£/L) 2 ) and using (£) = —ttc/QL 2 . However Eq. (32) is more 
general, and applies, for example, to the case of an interval in the interior of an finite 
open system (as considered in Ref. [IT]), or to a finite system at finite temperature. In 
fact at finite temperature we see that this correction (which is valid for £ <C /3) already 
exhibits the beginning of the crossover to thermodynamic entropy which actually 
occurs only in the limit £ ^> (3 (see Ref. [2] ) . 

Another case where 1-point functions can arise is in the case when the interval 
is near a boundary (say, at a distance ~ y) with f « 5. In that case the leading 
correction comes from the case where two of the {kj} are non-zero, equal to (say) k, 
and it involves the square of the 1-point function (<pk(y)) = fk y~ Ak ~ Ak , where is 
universal (but depends on the boundary conditions). The correction to Trp^ is then 

E^ 2, 4 ■ (34) 

3.2. Two intervals in an infinite system. 

We now apply the general result to the case of interest for the rest of this paper, i.e. we 
take A = \u\, vi] U {112, Vq\. Global conformal invariance gives Trp^ as in Eq. ([3]) with 
the harmonic ratio in Eq. Q. In terms of the notation above, we have l\ = \u\ — v\\ 
and £2 — \u2 — V2\. Let r + v\ — 1*2 — V2\ be the distance between the centers 

of the intervals. We are interested in the limit where ^1,^2 ^ r - I n that case 



l + O 



((W 2 )] , (35) 
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and the prefactor in Eq. ^ is 

7 ^~ U f>-^ = = ^(l + x + Otf)) . (36) 

{Vl - Ul){V 2 - U 2 )(V 2 - Ul){U 2 - Vi) Zit 2 {l-X) Zit 2 

Notice that often r is defined as the distance between the two intervals \u 2 — but 
the current choice absorbs some of the subleading corrections as clear from the fact 
that Eq. ( 35 1 has no corrections at order 1/r that would be present with any other 
choice of r. Furthermore if £% = £ 2 , we have x = £ 2 /r 2 exactly. 
From Eq. (16), we have 

■&pj = (/i/a)c», (37) 



and inserting the short interval expansion (17) for each interval, we find 

n 

tt pa = e e ^(WM^n^M^K^ 

{k l3 }{k 2j } 3=1 
n 

= E C {M(W^2)n<Mr)^0))c, (38) 
{**} i=i 
where in the second line we used orthogonality, i.e. the only non-zero contributions 
are from {kj} = {fei^} = {fc2j}[i|| Using now Eqs. (l22l and (l25| we have 



Trp\ = c^(44)- c/6(n - 1/n) E( £ i^) E ^ (Aj+A3) 4 fcj} n<^( r )^(0)>o 

{kj} 3=1 



fee \E,-( A i+ A ^) 



= 4(^ 2 )- c ^- 1 ^E(S&) <n<M e27rij/ ")>c- (39) 

Thus we obtain an expansion of Trp^| in powers of £i£ 2 /r 2 , which, via Eq. (35), can 
be converted to the expansion of the function F n {x) in powers of x. 

As explained above, to organize the sums as an asymptotic expansion (i.e. with 
increasing powers of x) we have to order all possible n-tuples of quasi-primaries fa. 
with non vanishing correlation functions according to the value of $^,-(Aj + Aj')- 

3.2.1. The contribution of the stress-energy tensor. There is an important term when 
one of the fa . = T (or T) and the rest of the kj = 0. Using the above values of dx 
and the fact that (T(r)T(O)) = c/(2r 4 ), this gives a contribution 

cn(l-l/n 2 ) 2 fh£ 2 \ 2 _ cn(l - l/n 2 ) 2 2 

72 ^ r 2 J ~ 36 X ' 1 Uj 

where the factor 2 in front of the lhs takes into account that there are two identical 
contributions from T and T correlations. Thus, in the infinite line, there are no 
contributions at order 0{x) in Trp^ coming from one-point functions. There could 
also be terms 0(i 2 ) arising from two fields each of dimension 1. 

|| This assumes that the whole system is critical and described by the CFT. In fact the short interval 
expansion can also be used in the non-critical theory, as long as l v -C the correlation length f, in 
which case these correlators arc proportional to the Zamolodchikov metric. 
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3.2.2. The leading contribution. The leading term, assuming that there are operators 
of sufficiently low dimension, comes from the case where two of the {kj} are non-zero, 
say equal to k. Their contribution in !F n {x) is 

s 2(A fe+ A fc ) J- (di JlJ2) f, (41) 
3i<h 



where d^ 1 is given by Eq. (29). There is an important point to be made here 



concerning the role of conformal spin Sfe = A& — A&. The expression for d^ 1 is 
periodic under (ji,j2) {ji + n,j2 + n), and therefore we can shift the summation 



variables (ji, J2) — > (Ji + I.J2 + 1) in Eq. (41 1 without changing the result. However 



^0'l + lj2 + l)^2 _ e -8^i(A fc -A fc )/n ^(jlh)^2 ^) 



Hence the sum in Eq. (41) vanishes unless Ask/n is an integer. This means that the 
contribution of a given primary field with non-zero spin is non vanishing only for a 
finite number of values of n (which are multiples of 4). 

From above, we can easily write down the coefficient Skin) entering in the 
expansion of the function T n (x) in power of x as 

. . p 47ri(ji+j 2 )sfc/™ 

s k (n) = 2 ix » Y, 

0<ji<h<n- 



\ e 2nij 1 /n _ e 2nij 2 /n\4x k 
i-l 1 1 



~~ ^ l^nMh - 3i)/n)\ 4xk ' 

where we defined Xk = A& + A&. Since 4sfe/n is an integer, the numerator is at most 
a sign. For Sk = 0, using the translational invariance in jx,2i the sum can be rewritten 

as 

n— 1 i 
n y—^ 1 

= 2^ (sinirj/n)^' (44) 

that reproduces the result in [18 . In any specific CFT, this contribution should 
be multiplied by the number of inequivalent correlation functions giving identical 
contributions (due e.g. to degenerate operators). 

3.2.3. The contribution of the three-point functions of primary operators. In the case 
of a three point function, the contribution in !F n {x) is given by 

a .A fcl +A fcl +A fe2 +A fc2 +A fc 3+A fc3 (4 jl fc 2 fc 3) ) 2 ' ( 45 ) 

0<ji<j2<j3<n-l 



where d^l^J is given by Eq. (|30|). Analogously to before, the contribution of non- 



zero spin operators vanishes unless 2(sk t + Sk 2 + Sk s )/n is an integer. We recall 
that in the following equations, to have a non vanishing contribution, we also need 
the coefficient c 12 3 to be non-zero (for example for the Ising model, the three point 
functions of the spin and energy operator are both vanishing, and the one with lower 
scaling dimension would be (a(zi)a(z2)e(z3))) . Thus, in most of the physical relevant 
cases, the contributions from three points functions are subleading compared to four- 
point ones, as we shall see in specific examples. 
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Specializing to the case of zero conformal spin, we have the coefficient entering in 
the expansion of the function !F n {x) in powers of x s^iji) = s klk2k3 (n): 

2 2(x 1 +x 2 +x 3 ) Cl23 



h<h<h l e " e " I 12,3 l e " 



E 



> — S — \2X13\2 \p — „ — 



Cl23 



(46) 



., ; ,v i| <„_Jsin(^ 21 /n)] 2 -i^3[ s i n ( 7rj31 / n )]2x 13;2 [ s i n ( 7rj32 / n )]2 a;2 3 ;1 

where j P q = j p — j g and x a b- c — x a + Xb — x c . Using the translational invariance in the 
jp, one of the three sums can be removed, but doing so is not very illuminating. 



3.2.4- The contribution of the four-point functions of primary operators. For a 
general CFT the four-point function has the form 

{<j) kl (z 1 )(j) k2 (z 2 )(j) k3 {z ? ,)(l) kl {z A ))c= ]J ( z p~ z q )^{z p - z q ) A ""G(x z ,x z ) , (47) 

l<g<p<4 

with A pq = Ylf—i Af/3 — A p — A q and similarly for A pq . x z stands for the four point 
ratio of the points z p . The universal function G(x z , x z ) should be calculated on a case 
by case basis because it depends on the CFT and on the specific operators considered 
[23 [23]. The function G(x z ,x z ) can be rather complicated for a general CFT, but in 
the two cases we will consider in the following and for the operators of lower dimension, 
it assumes a rather simple form that allows to write down d{ kj y explicitly. 

The four spin correlation function in the Ising model can be written as [231 123] 

1/2 r 



{a{z 1 )a(z 2 )a{z 3 )a(z 4 )} c = 



1 



2l 3 2 24 



Zl±Z23Zi2Z 34 



1 



212^34 



2l 3 2 2 4 



Z14Z23 



Z13224 



(48) 



that gives the lowest power in x coming from a four-point function in the short length 
expansion. For the free compactified boson, the lowest contribution is given by the 
correlation function of four vertex operator two with charge y/JL (V^z) = e 1 ^ 1 ^ 2 )) 
and two —^ffL {V-^(z) = e~ l ^ v ( z )). The value of fj, = a = min[?7, l/rj\ is directly 
related to the compactification radius. The three relevant correlations are 

2p 

<V^i)V^ 2 )V_ M (23)V%(*4))c = 



212 Z34 



244223^14223 

2l32 2 4 



214223212234 
223214 



2// 



2/j 



(49) 



<M^-^2)V^(2 3 )^(24)>C = 

224213212234 

The correlators with fi — > —\i should be also included and these give identical 
contributions. The correlations for the Ising model and for the free boson can be 
treated on the same foot since the formers come from the bosonization of the latter 
with \x — 1/4. 

It is now easy to obtain the factors d^ll\^^\ with [i p = ±/i as for example 



(V p {e 
(2n)- 4 ^ 



sin(7r j 2 i/n) sin(7r j 43 /n) 



•)>< 



2/' 



sin(7rj42/n)sin(7rj3i/n)sin(7rj4i/n)sin(7rj 32 /n) ' ^° 
and the other two just by permuting properly j p . We used that all these operators 
have zero conformal spin. 
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4. The Ising model 

The partition function of the Ising model on a generic Riemann surface of genus g 
with period matrix r has been obtained in Ref. [35] (see also and reads 



Z=(2 qu ) 1 /2 2 -sV 



e 



(0|r) 



(51) 



where Zq U is the quantum part of the partition function for a free boson that does 
not depend on the details of the Riemann surface (and on the target space as well) 
but only on its topology. Zq U has been calculated by Dixon et al. [27] (see also [28 J. 
In Eq. (51), the sum is intended over all the 2 9 vectors e and 6 with entries and 



1/2. We mention that this result for the Ising model is rather special, indeed usually 
it is not possible to give a closed formula for a general Riemann surface with arbitrary 
period matrix t for all CFTs. For example, also for the compactified boson, such a 
formula does not exist (in [I], the property that the matrix t is purely imaginary has 
been exploited to obtain Eq. (JHJ). 

Combining Eq. (51) with our results in Ref. [T], it is easy to derive the Ising 
partition function on the n-sheeted Riemann surface TZ n ,2- In fact, a by-product of 
the calculation for the free boson pQ is that the (n— 1) X (n— 1) period matrix is given 
by Eq. pi). Although derived for a free boson, the period matrix is a pure geometrical 
object and it is only related to the structure of the world-sheet TZ n .2 and so it is the 
same for any theory. 



The function F n {x) for the Ising model is proportional to Eq. (51) with r equal 
to the matrix T in Eq. Q, and the proportionality constant is fixed by requiring 
•F n (0) = 1- In this way, we simply obtain Eq. (11). We mention that for known 
formulas analogous to Eq. (51) and valid for other theories, T n {x) can be easily 
obtained in the same way. 

We performed several checks for the correctness of this equation. For n = 2 it 
reduces to 



1 

71 



1/2 



+ x 1/4 + [{l~x)x} 1/4 + (l-x) 



1/4 



a/4 



(52) 



obtained (and checked against numerical calculations) in Ref. [15]. For n = 
3, 4 numerical estimates of these functions have been obtained from the exact 
diagonalization of the Ising chain in [17) . The comparison of our prediction with 
the numerical data for n = 3 and n = 4 is reported in Fig. [T] The agreement is 
excellent, taking also into account that the extrapolation to £ — > oo in Ref. [T7] has 
been obtained without knowing the asymptotic result. 



4-1. Invariance under x — > 1 



When x — > 1 — x, T n (1 — x) is given by Eq. ( 11 ) where T is replaced by T given by Eq. 
^ where x is replaced by 1 — x (i.e. (3 y — > l//3 y ). We denote by J- n { x ) — J~ n {l — x) the 
expression obtained by using Y instead of Y in (11). In order to show the invariance 



of ( 11 1 under x — > 1 



e 



x, we observe that in general we have 

„2m a t -(z+b) 



y/det(-iT) 



9 



z + b 
—a 



(Ol-T- 1 ), 



(53) 
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Figure 1. ^-3,4(2;): The extrapolated numerical data of Ref. I17| are compared 
with our prediction. Only x < 0.5 are reported because of the symmetry x — > l—x. 



for all complex vectors a,b and z. This can be shown by a straightforward 



generalization of the first part of the appendix B of Ref. [T]. The relation (53) 
implies that 



1 



6(0|T) 



9 



a 
b 



(0\T) 



1 



9(0| -T- 



e 



(0| - T~ 



(54) 



Applying this relation to T = T = —AT 1 A (the transformation matrix A is defined 
in App. B of Ref. pQ), we find 



F n {x) = 



2 n - 1 @(0\A- 1 TA) 



(OlA^TA) 



(55) 



In Ref. pQ we have used that Q^A^TA) = 6(0|r) without proof (see Eq. (B.5) of 
PQ). Here we argue that a generalized version of this equation is 



e,6 



^2\e d £ (o\a-'ta)\=^2\q 



e,S 



(o|r) 



(56) 



whi ch is numerically true and we give without proof. From this relation, we find that 
(11 ) is invariant under x 1 — x, namely J- n (l — x) = J- n {x) = F n (x). 



5. The small x expansion 

As we showed above, for small x, the scaling function J-"„(a;) displays a systematic 
expansion revealing the full operator content of the theory in an explicit manner. 
Here we consider the first two non vanishing contributions from the expansion of 
T n {x) for the free boson and the Ising model. This expansion has the form reported 
in Eq. (10) that gives the first terms in the asymptotic expansion as long as 2a < 1, 



because 0(x) terms exist as well (the latter are actually trivial as we shall see and do 
not correspond to any operator). We recall that a is the double of the leading scaling 
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operator in the theory (i.e. for the boson a = minfry, l/rj\ and for Ising a — 1/4). 
The coefficients s% 4,(11) are explicitly calculated in the following. We mention that in 
both models the contribution of non vanishing three-point functions are subleading, 
as clear from the conformal spectra of the two theories. 

Notice that, since the exponents in this expansion do not depend on n, the 
structure of the power series passes over to the analytic continuation and so the scaling 
function of the von Neumann entanglement entropy is 



VN\ 



(*) = = Q a 4(1) + (|) 2Q 4(1) + • • • , (57) 

where we used .Fi(x) = 1 (i.e. TrpA — 1). The calculation of the analytic continuation 
of the amplitude s m (n) and their derivative at n = 1 is rather complicated and will 
be considered in Sec. [H 

5.1. First order for the compactified boson 

The calculation of the leading coefficient S2 (n) for a free compactified boson has been 
already reported in Ref. [T]. This is repeated here for completeness and to fix the 
notation for the Ising model and for calculating the second non vanishing order. 
The first ingredient we need is the expansion of Pk/n f° r 00 ~* 

p k = _ s ^±R( logx + fk+ y f k ^2 lE + ^( k -\+^(i- k 




n J \ n 

where is the Euler gamma, ip(z) = ^'(z)/T(z) is t ne polygamma function and 
Pi(z) is a polynomial of degree 21, whose explicit expression is not needed because it 
contributes at order o(x) that we do not consider here. Plugging this expansion in the 
Riemann-Siegel theta function we obtain 

0(O|7yr) = 1 + x^e^™ [1 + O(x)] , (58) 

m£Z"- 1 \{0} 



where the two vector-dependent coefficients a m and b m are 

^ sin 2 7T- jm^Cfc/n-m, (59) 



n 
fe=i 



E sin2 



n 
fe=i 



t- ) /fc/ntti • C k/n m, (60) 
and the matrix Ck/ n has elements pQ 

( C 



27T — (r — s) 
n 



(61) 



The leading term in the small x expansion is provided by those vectors m £ Z™ 1 \ {0} 
which minimize a m . When evaluated on the vectors m, it is easy to show that 

n— 1 n— 2 ^2 1 ^2 ^ 

im = ^m j 2 -^ m i+1 = 1 2 + - ^(toj - m J+1 ) 2 , (62) 

3=1 j=i i =1 

where indicates the j component of the vector m. When minimizing a m , the 
second term in the rhs of the last expression tends to make all the components of m 
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equal and only "domain walls" give a non vanishing contribution, proportional to the 
square of the difference of the two neighbor components at the wall. The first term is 
like a "boundary term" and tend to make the absolute value of of the boundaries mi 
and m„_i as small as possible. Within these ingredients, it is obvious that (among 
the non zero vectors) a m is minimized by all vectors of the form 

;0,...,0,±1,...,±1,0,...,0), (63) 



4 



h h-h 

with < jx < 32 < n — 1 that give all a m = 1. In order to calculate b 

' sin (7r%(j2 - jx)) 
' pq \ sin (ir 



/ \ 9 

32-31 ' ' • ~ 



s(h - h) - ml • C k/n ■ m ± = E N = ^r^V ^ ■ ( M > 

^— ' V "/ pq \ Sin (7T- I / 

and the integral representation of the polygamma function 

^(2/) + 7£= / -r-dt. (65) 

Jo i — e 

Exchanging the order of sum and integral and then performing the latter, after simple 
algebra we find 



b m± = -2 log 



,32- Ji 
An sin 7r 



(66) 



We now can use the translational invariance of s(j 2 — jx) to fix jx = and multiply 
by the number of possible jx at fixed j2 — ji, i-e. (n — j-i)- There is another factor 2 
to take into account the two possible signs in the vectors m± . 

Now, assuming r\ ^ 1, in the ratio of Riemann-Siegel theta functions occurring in 
T n [x) the leading term is given by the minimum between r\ and 1/rj, that we indicate 
as a — min(?7, 1/rf). Therefore we get the small x behavior of the scaling function 
(renaming j 2 = j) 

■ = i [2nsin(7ri)J V4n 3 =i [sin (tt^)J 



n— 1 , 

,-=i [sm(7ri)J 

where the dots denote higher order terms in x. We stress that this expansion is valid 
only for r) ^ 1, in fact for n = 1 the denominator in J- n (x) (that is of order O(x)) 
cancels exactly the numerator. 



This first order term in Eq. (67) coincides exactly with the prediction from 
short length expansion ( |44[ ) with Xk = a/2 corresponding to the scaling dimension of 
the vertex operators V± a (z). The additional factor 2 comes from the fact that both 
correlations (V a V- a ) and (V- a V a ) must be considered. 

5.2. Digression on the denominator 

The previous analysis also provides the small x expansion of the denominator in T n (x) 
corresponding to a = 1. We have 

n— 1 1 

eWW + n^V „ +-. (68) 

3=1 l s ™{iri)] 
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As showed in Ref. [T] we have 

n-1 



1=1 



1 1, 9 
2 = o( n - !) 

sm(7r?/n)] 3 



so that the denominator is 



[0(o|r)] 



G 



0(x 2 ), 



(69) 



(70) 



that gives an 0{x) contribution to the scaling function !F n {x). However there is no 
operator in the spectrum of the free boson whose correlation function can contribute 
at order 0(x). This fact could sound as a breakdown of the short-length expansion, 
but it is not. Indeed, it predicts the asymptotic expansion of the full Trp^ and in 
general any term in this expansion is also reproduced in J- n (x). However, in the case 
of integer powers, we should properly take into account the prefactor of J- n {x) in Trp^ 
as in Eq. ([3]). According to Eq. (36), this is (for c = 1) 

-| l/6(n-l/n) 

(71) 



1 



l + x-(n-l/n) + 0{x 2 ) 



(^ 2 )V6(n-i/n) 

that exactly cancels the 0(x) contribution of the denominator in Eq. ( 70 ). Thus there 
is no 0(x) contribution in Trp^ in agreement with the short-length expansion. 



5.3. First order in x for the Ising model 

We can now move to the calculation of the first non vanishing contribution to Tnix) 
for the Ising model. In order to study the behavior of Eq. ( 11 ) when x — > we isolate 
the contribution of the terms with e = 

1 u (o|r)|+££|e 

E^O S 



F n (x) 



2 n - 1 e(o|r) 

l + 0(x 



o 

6 



(0\T) 



)n—l 



£[i + o(*)]+£E| 

S e^O S 



(o|r) 



(72) 
(73) 



By using the definition of a m and b m in the previous subsection, we find for e =/= 



e 



(0|r) = E e 27rl(m+£)t ' s x a ™+'e h ™+' [1 + 0[x)] 



(74) 



m£Z" 



Again the leading contribution is given by all the vectors m that minimize a m . By 
simple inspection, the minimum value of a m + e is 1 /4 and for any e of the form 

£' = (0,... ,0,1/2,... ,1/2,0,..., 0), (75) 



31 32- jl 

the value a m+e = 1/4 is obtained for the two vectors 




,0) 



m 



(0,...,0) 



(m 



(m 



(76) 



(77) 



where m± are the vectors employed for the free boson in Eq. (63). It is therefore 
clear that we can exploit the results of the free boson to obtain the expansion for Ising 
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very easily. In particular, a m + e and 6 m + e are 1/4 of the corresponding free boson 
quantities. Thus from Eq. (l73|) we have 



T n {x) 



1 



2 n-l 
= 1 + - 

= 1 + 



^EEl 

£^0 S 

I — l" E E E 



(o|r) 



e 27ri(m+e) t -(5 :z ,a 



<5 m£Z" 



,1/4 



2"- 



E 



= i 



= i 



0<ii<j2<n— 1 
n-1 

3=1 



X 



1/4 



8 



+ 



,1/4 



n-1 

2n-i ^ ^ 



( n - 3) \ - 
2n sin(7rj/n)J 4 s 



(78) 
(79) 
(80) 
(81) 
(82) 



The various steps are very easy, but it is pedagogical to summarize the properties we 

_£ l± . is real and independent of 8, ± and j\. In (811 

In 



employed. We used that 6( m+ . 

we have used the translational invariance to fix j\ = and renamed j2 = j 
we have employed the free-boson result exposing the role of the factor 1/4 which is 
relevant in this case. The next step is to observe that • S = q/2, where q is the 
number of 1/2's occurring in the first j coordinates of S. The remaining n — 1 — j 
coordinates of d do not change the scalar product and this freedom provides a factor 
2«-l-J. Thus we find 



S 

The final result is 



F n (x) = 1 



= 2 n - 1 - j 



O - i) 



3 

E 

9=0 



,TTiq/2 



+ e 



—iriq/2 



= 2 % 



;i E 

j=\ 1 2n sin(7rj'/n) 



• = 1 



— 1 
4n 2 J 



E 



sin 2 [7rj jn\ 



(83) 



(84) 



Notice that this equation has the same form as the one for free boson Eq. (67 1 for 
a = 1/4, with an additional factor 1/2 in front. This is exactly what predicted by 



short length expansion (cf. Eq. (44)) from the two-point spin correlation function. 



As for the free boson, the 0{x) contribution from the denominator of T n (x) 
cancels with the 0{x) contribution of the prefactor in Eq. (|3j), in agreement with 
the short length expansion, since also for the Ising model no correlation function can 
contribute at O(x). 



5.4- Second order for the compactified boson 

To obtain the second order in the small x expansion, we need to know the vectors 



m in Eq. (58) that give the second minimal value of a m . Since a m can assume only 
integer values (cf. Eq. (62)), we need to identify all vectors giving a m = 2. From the 
rhs of Eq. (62 ), we can have a m = 2 only by introducing another two domain- walls of 
height one in the string of ±1 in m± in Eq. (63). Thus the vector m giving a m = 2 
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can only be of the forms 



0. 



,0,±1,...,±1,0,...,0,±1,...,±1,0,...,0), 



(85) 



3l 32-Jl 33— 32 34-33 

m\ EE ( 0, . , Q, ±1, . , ±1 ,0, . , 0, T l, . , T l , 0, . . . , 0) , 

31 32—31 33—3*2 34—33 

n4 = (0, . 0, ±1, . ,±1 , ±2, ,±2 , ±1, . . ,±1 ,0, . . . ,0) , 

3i 32-31 33—32 34-33 

where < ji < j% < j 3 < j'4 < n — 1. For each m r we have two types of vectors 
corresponding to the sign of the first non-zero element (when needed we will specify 
the vector with an additional ± subscript). We remark that these vectors exist only 
for n > 4. For n — 2 and n — 3 the minimum value of a m greater than 1 is a m = 4 
and a m — 3 respectively, as predicted by short-length expansion. 

In order to get b m for the vectors m r above, the first step is to compute 



m • C, 



k jn 



m r . We introduce the sum of the elements of the matrix CV/ n enclosed in 
a rectangle made by A r rows and A c columns starting at the (p+ l)-th row and at the 
(q + l)-th column. It is given by 



p+X r q+\ c 



S(P, K,Q, A c ) = ^ ( C k/n)i3 
i=p+l j=q+l 

sin(7r(/c/n)A r ) sin(7r(fc/n)A c ) 



sin(7r(fc/n)) 2 



cos ( 7T— (2[p — q] + X r — X c 
n 



(86) 
(87) 



1 



2 sin (irk/n) 



sm 



(p-q+ X r ) 



sm 



nk 



(jp-q- X c 



nk 



(J) - q + X r - A c ) 



nk 



(p-q) 



In the special case of a square (i.e. X r — X c = A) on the diagonal (i.e. p = q) it reduces 
to Eq. (64), i.e. s(p, X,p, X) = s(A). Then, we have 

m' • C k j n ■ m = s(j 2 i) + s(j 43 ) + 2s(ji, j 2 i, j 3 , ks) , (89) 

= s{j 2 i) + s(j 43 ) - 2s(ji, j'21, 33,343) , 

m 2 = s(j 21 ) + s(j 43 ) + 2s(j 3 , J43, ji, 321) 

+ 4(s(j 32 ) + s(j 2 ,j 3 2, 33,343) + s{ji,j 2 i,j 2 ,j 32 )) , 

where we introduced jki — jk — ji ■ 

At this point we have all the ingredients to calculate b mr from their definition 



m\ -C k/n - mi 



m* • C k/ 



(60 1 and from Eq. (65) to deal with the polygamma function. It is a tedious, but 



straightforward calculation to obtain these coefficients. After long algebra, we finally 
have 



log 



4n 2 



where we defined 



for r = 0,1,2. 



sin(7rj 42 /n) sin(7rj 3 i/n) 



Qi 



sin(7rj2i/n) sin(7rj 43 /n) sin(7rj 4 i/n) sin(7rj 32 /n) ' 

sin(7rj 4 i/n) sin(7rj 32 /n) 

sin(7rj 21 /n) sin(7rj 43 /n) sin(7rj 42 /n) sva{nj 31 /n) ' 



(90) 

(91) 
(92) 
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g 2 _ sin(7rj 2 i/n) sin(7rj 43 /n) . 

sin(7rj42/n)sin(7rj3i/n)sin(7rj4i/n)sin(7rj 3 2/n) ' 

We have Q r > for any r. Notice that <5o,i,2 are invariant under translation 
(i.e. depend only on the differences of the j's) although Eqs. ( [89] ) are not. There 
are few other important symmetries of these factors. Each Q r is invariant under 
1 H 2 and 3 O 4. Under the cyclic permutation {1,2,3,4} O {4,1,2,3}, we have 
Qo - > — Qo an d Q1.2 — > —Q2.1, but the minus signs do affect the value of 6 mj . in Eq. 



(90 1 since they are function of the squares of Q r . 

After having computed the amplitudes b mr from the vectors m r , we need to sum 
over all of them according to 

s 4 (n)=2 ]T [Ql a + Ql a + Ql a ] , (94) 

0<Jl<j2<33<J4<n-l 



where we factorized out the factor (An 2 )~ 2a according to the definition (10 1 and we 
used the positivity of the factors Q r (otherwise the absolute value of each Q r should 
be used). The factor 2 comes from the two possible vectors (i.e. starting with ±1) of 
each type. 

This result for 84(71) agrees with the prediction of short-length expansion, using 



Eq. (50 1 and summing over j p . The three Q r correspond to the three possible four- 
point correlation functions of the vertex operator V± a (z), as clear by comparison. The 
factor two in front corresponds to the two sets of non vanishing correlators. 

It is also possible to further simplify 84(11). We can use the invariance under 
translation of the factor Q r . This implies that if we, for example fix j'2,3,4, the sum 
does not depend on ji that we can fix to and we get an overall factor a factor n — j'4, 
leading to 

s 4 (n)= J2 2(n-j4)[Q 2 a + Qi a + Qt a ] (95) 

l<j2<j3<J4<n-l 

= \ E [Ql a + Ql a + Ql a \ = l E [Qo a + 2Qn (96) 

1<J2 <fa <ji<n— 1 1<J2 <J3 <ji<n— 1 

In the last line we have used the identity 

E Qi a = E ^2", for any a, (97) 

l<j2<33<34<n — 1 l<ia<j8<J4<T» — 1 

that is valid only at the level of the sum (and after fixing j\ = 0), but not for the 
single term. This is a consequence of the permutation symmetry of the Q r 's. Another 
interesting relation, that is particularly useful in the case a = 1/4 (and so for the Ising 
model in the following) is 

QT = Q\' 2 + QT , (98) 

which can be easily shown by using the prosthaphaeresis formula involving the product 
sin A sin B. 

5. 5. Second order for the Ising model 

We can employ the results of the previous subsection to compute the second non 
vanishing term in the expansion of J^ n (x) for the Ising model. For a given characteristic 



Entanglement entropy of two disjoint intervals in CFT II 



20 



(e, 8), the vectors m £ Z" 1 contributing to the 0(x 1 ^ 2 ) order are characterized by 
a m+e = 1/2- From Eq. ([79]), we have 



s 4 (n) 



(4n 2 ^/2 



g 2'7ri(m+£) t S gbm+e 



(99) 



-EE | E 

e#0 5 ms.t.a m+c =l/2 

The vectors m contributing to this sum (i.e. with a m + £ = 1/2) are such that m + e 
are half of the corresponding ones for the free boson, i.e. half of m r introduced in Eqs. 
(85 1. Thus, as for the boson, the sum over m can be rewritten as sums over jx,2,3,4) 
but the choice of m also fixes e, and so this sum can be dropped. Using the value of 
b mr in Eq. ([90]) we have 



E 



— y 

n-l 



(100) 



0<Jl<j2<j3<j4<n-l 



Qq 2 cos(7rm • 8) + Q{ /z cos(7rnii ■ 8) + Q 2 ^ jz cos(7rm 2 ■ 8) 



a/2 



V2, 



To write the sum over 8 in a manipulable way, we notice that the scalar products are 
very simple to write as 



m • 8 = 



qi + <?2 



mi • 8 = 



<7i - 92 



m 2 • 8 = 



qi+q 2 + 2q 



(101) 



2 2 ' 2 

where q\ is the number of 1/2 in 8 in the part of m r , q 2 the number of 1/2 in 8 in 
the part j'43 of m r , and qo the number of 1/2 in 8 in the part j'32. The argument of the 
sum does not depend on the presence of 1/2 in the j'-th component of 8 with j > j\ 
and j < ji, so that the 2 n ~ 1 ~^ il possible choices give equal contributions providing 
an overall factor 2™ -1 ~ : ' 41 in front of the sum over (ft), 1,2- Putting everything together, 
we have 

J21 J43 J32 

s 4(«) = E E E E 

0<jl<j2<j3<j4<n-l 81=0 92=0 qo=0 

x cos [|(gi +<7 2 )]Qo /2 +cos [-(gi - g 2 )] Qj /2 + cos [|. 

Since, for any triplet of integers {gi, q 2 , qo}, the parity of qi+q 2 , Qi — q2 = (<7i+<72)-2g 2 
and qi + q 2 + 2qo is the same, the terms having odd qi + q 2 do not contribute. Thus, 
we have 

321 343 332 



321 

qi 

1/2 



J43 

<?2 
COS 



.732 

'7o 



(102) 

( qi +q 2 + 2qo)]Ql /2 



s 4 {n) 



JZl jq^-i JO^ 

^ ' 2/41 — 1 ^ < ^ < ^ < 



0<jl<j2<j3<j4<n-l 



91=0 92=0 9o=0 



./2 1 
'71 



.743 
<?2 



.732 

'7o 



C0S UWl + 92. 



g /2 + (-i) 92 Q; / ' + (-i) 9o g 



1/2 



,1/2 



(103) 



(104) 



where | cos(-7r(qi + q 2 )/2)\ simply selects the terms with even q± + q 2 . Thus, the sums 
over q's in (|103[) can be organized as follows 



EEE 

91 92 90 



even 91+92 



E +E 

even 91 odd 91 



E +E 

even q 2 odd q 2 



even q\-\-q-2 



E E +E E 

even 91 even 92 odd 91 odd 92 " 

E + E + E 



E + E 

even 90 odd 90 

+ E 



E ( 105 ) 

90 

(106) 
(107) 



91,92, 9o£2Z 9!,g 2 e2Z gi ,g 2 e2Z+l 91 ,92 ,9o£2Z+l 
9o£2Z+l 9o£2Z 
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In Eq. ( 105 1 the condition of even q\ + qi selects two of the four sums coming from 



■iqi 



^2 J2 Q2 - Each of the four sums in (107) is characterized by a specific sequence of 



relative signs between the three contributions to the general term in (104 1. Once the 



relative signs have been fixed, the dependence on g's is only in the binomials. This 
means that we can perform the three sums over g's in each one of the four terms 



indicated in (107). In order to make this sum, we find it convenient to leave the parity 
of three indices g's unrestricted and introduce a factor | cos(7T(7/2)| or | cos(7r(q + l)/2)| 



when q is required to be respectively even or odd. Then, we have that 

j 



3 

E 

9=0 



cos 



E 

g=0 



cos 



r(? + l) 



(108) 



By using (107) and (108) we can perform the sums over g's in (103) and reduce them 



.1/2 



to the sum of four terms. Thus, we have 

«4(n)= £ jE|Qo /2 +^; 

0<Jl<j2<j3<i4<n-l n,T 2 

where ti,t 2 € {—!,+!}. Now, we use Eq. (|98| which allows us to further simplify 



r 2 Q 2 



(109) 



the sum over (ri, T2) in ( |109[ ), which becomes 

E 



il' 2 



T!Q\ /2 +T 2 Q 



1/2 
2 



= 2 



(\qI 



1/2 



f/ 2 



1/2 



(110) 



1 /2 

where we have also used that Q ' > 0. Notice that, because of Eq. (98), the term 



corresponding to (ti,T2) 
simple as 

s 4 (n) = 



( 



E 

0<ji<i2<j3<i4<n 



) vanishes identically. Thus, the final result is as 

(111) 



.1/2 
so 



that is one quarter of the free-boson result for a = 1/4. This value of Si(n) agrees with 
the result from short-length expansion using the four-point spin correlation function 



in Eq. ( 48 ) and evaluating it at the roots of unit as in Eq. ( 50 1 . This would give the 



sum of the three pieces in Eq. (110), that we have shown how to further simplify. 



As for the free boson, we use translational invariance to fix j\ — at the price of 
the overall factor n — j^, to rewrite s^in) as 

)l/2 _ 2 



S4,(n) 



l<j2<j3<J4<n- 



(n - j 4 )Qo 



E 

I<j2<j3<j4<n- 



Ql /2 - 



(112) 



6. The analytic continuation 

In this section we consider the analytic continuation of S2(n), and in particular 
we obtain its derivative at n = 1 giving the small x behavior of the entanglement 
entropy. In Ref. pQ we provided a general method allowing a very precise numerical 
computation of S2 (n) for non-integer n, but we have been able to give a full analytic 
answer only for specific values of a, i.e. for a = 0, 1/2, 1. This has been generalized in 
Ref. [IB] to any integer multiple of a/2. Here we fill this gap and give a full analytic 
answer for the analytic continuation for any value of a. 

We start rewriting S2 (n) as (we temporarily drop the factor Af) 

n-l 



S2(n) 



1 



(113) 
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and since it is periodic, let us write a Fourier series 

(sin 2 (u/2))- a = ^/(fc)e 4fc ", (114) 
kez 

where 

f(k) = — [ 27r {sm 2 (u/2))- a e~ lku du = - Asm 2 y)- a e - 2iky dy . (115) 
2tt Jo n Jo 

For a < 1/2, this integral can be done (see e.g 3.631 of Gradshteyn and Ryzhik f2~§\ ) 
, . 2 2a e- l7Tk T(l - 2a) , s 

f W = ra- a + k)r\i- a -ky (116) 

Note that f(k) ~ |fc|~ 1+2a as \k\ — > oo, as expected by power-counting. 
We need now to perform the sum 

n-1 



which gives n — 1 if k = (mod n), and —1 otherwise. The sum for S2{n) can be then 
written as 

«a(n) = i -/(A;)], ( 118 ) 

kez 

that vanishes at n = 1 as it must. This formula could provide the required analytic 
continuation, but it should be handle with a lot of care. Indeed there are two main 
problems with this sum, one numerical and the other conceptual. 

Let us start with the first problem. The two subsums defining S2(n) (i.e. nf(nk) 
and 5Z/(fe)) are both divergent. Their sum can be made convergent only if the two 
terms are summed in the proper order, i.e. if we use the same upper limit for nk in 
the first and k in the second and only within this recipe we can take the difference and 
send the cutoff to infinity. If one makes the sum naively, it diverges. For integer n, 
by summing numerically in this way, we reproduce easily the correct results obtained 
by the direct sum Eq. (113). However, for non- integer n the obtained sum is not 
a smooth curve (see Fig. |2|, signaling some troubles with this analytic continuation 
that leads to the second problem. 

Indeed, it is obvious that the analytic continuation of a function defined only 
on integer numbers is not unique, because, e.g. we can multiply it by some function 
like cos 2 (7m) (that is 1 on the integer) or add any analytic function times sin(7rn) 
that vanishes for integer n. However, all these analytic continuations with the same 
values on integer numbers have different behavior at ±ioo and this is the essence of 
the Carlson's theorem [3D], stating (in physical words) that the analytic continuation 
is unique once the desired behavior at ±ioo is imposed. While at first sight this 
can seem a problem, it has been indeed shown that this arbitrariness can be a great 
advantage. For example, it can provide the solutions of different physical problems 
with a single result for integer numbers, as e.g. done in Rcf. |31j where correlation 
functions of one-dimensional models with different statistics have been obtained by 
the same replicated form, only requiring different behavior at ±ioo. 

Going back to our problem, it should be clear that to select the correct analytic 
continuation we should have a look to n — > 00, i.e. we should consider the behavior of 
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f(k) for k — > oo. Exploiting the asymptotic expansion of the T function [29j . we have 
f(k - ±oo) * 2-e--r(l - 2« ) Sil ^|M (119 ) 
= 2 2 «(cos(.fc) - ism(7rfc))r(l - 2a) M^)coB(7rfc) + 8m(|fc|7r)coB(7ra) 

In this simple formula we see in action all the pieces for which Carlson's theorem has 
been formulated: there are additive sin(7r/c) pieces (that do not change the sum for 
integer k) and multiplicative cos 2 irk (that just give 1 on integers). Both factors change 
the value of f(k) [and so of S2(n)] only for non-integer values. Carlson's theorem states 
that we need to impose the required analytic behavior for ±ioo to select the correct 
analytic continuation. In our case we require that S2(n) (and so any term in the sum 
f(nk)) does not blow up exponentially. This implies that all the sin(7rfc) terms are set 
to zero and the cos 2 (7rfc) to one. In order to do so, we perform the replacement 



/(/•- ±m> ,/(/.■ ±x ) = 2-'Tl I -2a) (120) 



f(k — > oo) and g(k — > oo) have the same value for integer k, but they arc different for 
non-integer ones. 

Being guided by the asymptotic behavior for large k, we can employ the same 
trick for the full f(k) in order to have the correct analytic continuation. We perform 
the replacement 

= /(fc)sin( TO )e^ = 2 2 "r(l~ 2a )sin( TO ) 
n > yK >~ sin[^(a + |fc|)] sm[7r(a+|fe|)]r(l-a + fc)r(l-a- k) K ' 



4 Q sin(7ra)r(l - 2a) sin[7r(a - \k\)]T(a - k)T(a + k) 



(122) 



Again, f(k) and g(k) have the same values for integer k, but it is easy to check that 
while f(k) oscillates for non integer k (signaling 'bad exponential' behavior at ±ioo), 
g(k) is monotonic. 

The sum for S2(n) can be rewritten as 

oo 1 oo 

fi 2W - 2 E [ng{nk) ~ g{k)] = ^H (0) + EN(«fc) - gm , (123) 

k— — oo k — 1 

that is the analytic continuation we were searching for, but valid only for a < 1/2. As 
stressed above, this sum is convergent if for the first and second piece that same cutoff 
is used and after sent to oo. We performed numerically this sum for non-integer n 
finding perfect agreement with the results in Ref. [1 obtained by Borel resummation. 
Fig. [2] shows the difference in the summation of s^in) obtained for non-integer n while 
using the function g(k) instead of f(k). 

Eq. (123) allows analytic manipulations that give the exact result for n — 1. Let 
us use the integral representation of the product of T functions 

T(p)T( q ) = 2T(p + q) rVs^-Vn^) 2 *- 1 , (124) 
Jo 

to rewrite 

= 4" sin(^)r(l - 2a) _ [fc m 2r(2a) /"^ ^2a-2fc-l (gill0) 2a+2fc-l # 

7T 2 Jo 

2 2Q sec(7ra) . /' 7r/2 / sin2(/) N 2Q-i , 1 x 2 fc,j ,,„ r , 

= -sin[7r(a- \k\)} / (— ^) (tan^) 2fc d</> , (125) 
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— a=0.4, g(k) 
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■ a=0.1,g(k) 











r. 1 I I I I I 

1 1.5 2 2.5 



Figure 2. This figure shows the importance of the Carlson's theorem to choose 
the correct analytic continuation. For two values of a (0.1 and 0.4), we report 
Sn/(n> — 1) as function of n, as obtained by summing with the function f(k) and 
g(k) (i.e. using Eq. ( | 1 1 8 [ ) or Eq. ( |123| l respectively). It is evident that when using 
f(k) an oscillatory behavior (with period 1) is superimposed to a smooth curve, 
while for g(k) the resulting S2(n) is smooth. For integer numbers > 2 they have 
the same value, as they should. The limits for n — > 1 are very different and the 
one obtained from g(k) is the correct one. 



where we used also 4 Q r(l - 2a) sin(a7r)2r(2a)/7r 2 = 2 2a sec(7ra)/7r. 

In order to evaluate y\ _-, ng(nk), we exchange the order of sum and integration. 
Calling v — tan <fi, we need the sum 

v 2n [sin(7r(o; — n)) — sin(7ra)w 2 ™] 



I n = sin[7r(o! — nk)]v 



2k n 



k=l 



v 4n - 2cos(nn)v 2n + 1 



(126) 



that is strictly valid only for |k| < 1 (i.e. <f> < 7r/4), but we will assume its validity 
everywhere. It is easy to take the derivative of nl n wrt n and set n = 1, obtaining 



d{nl n ) 



On 



[n cos(7ra) — (1 + v 2 + 2 log v) sin(7ra)] 



=i (l + v 2 ) 2 ' 

= sin 2 (0) cos 2 (4>) [7rcos(7ra) - sin(7ra) [sec 2 (j) + 2 log(tan <f>)] ] . (127) 
Plugging the various pieces together we hav^j] 



2 S ' 2 (1) = 5 (0) 



1 



7r cos(7ra) 



7T/2 



(sin20) 2a+1 7rcos(7ra) 



% The following integrals are useful for this calculation: 

rr/2 



(sin 2(f) 



2c+i = y^rca+i) 
2r(a+|) 



tt/2 



(sin 2cf>) 2a+1 (2 log tan <f> + sec 2 <j>) 



AFr(a) 4 a r(l - 2a) tt 



r a 



r(l — a) 2 tan ira 
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,ir/2 

- / (sin20) 2Q+1 sin( 7 ra)(sec 2 + 21og(tan< ? (>)) 
Jo 

_ 4 a r(l - 2a) y/^T(a + 1) tan(7ra)r(a) _ 0iT(a + 1) 

= r(i-a)2 + 2r (a + l) ~ V^r ( a + f ) = 2r (a + 1) ' (128) 

that is the result anticipated in Eq. ( [l4| , where we restored the multiplicative factor 
A/". While this result has been derived for a < 1/2, the final form of s 2 (l) is analytic 
all the way up to a — oo and indeed it agrees both with the numerical results of Ref. 
[T] for all a < 1 and it is straightforward to show that for positive integer a reproduces 
the right result (in this case, s 2 (n) is a polynomial [TB] and it is easy to find s' 2 (l)). 



6.1. Consequences for the Ising model 

In the case of the Ising model, the first order in x of the scaling function J- n (x) is half 
of the result for free boson at a = 1/4 (cf. Eq. (84)). Thus the scaling function for 
the entanglement entropy is 







'5\ 



Fvn{x) = ^n{x) t = 7 ^\\\ x 1/A + 0(a: 1/2 ) - 0.309012 . . . x 1 / 4 + • • • . (129) 



4V2r (I 



Note the numerical value of the amplitude of the term x 1 / 4 is close to (but definitively 
different than) 1/tt = 0.318 . . ., that was proposed in Ref. [15 on the basis of purely 
numerical results. We checked that the numerical precision of the data in Ref. [T5] 
does not allow to distinguish between these too close values. 
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Appendix A. Mapping of two intervals to the cylinder geometry 

Some of the discussion of the main part of the paper may be more readily 
comprehended if we conformally map the n-sheeted conifold consisting of n planes 
connected along the two intervals, to another consisting of n open cylinders connected 
in a certain way along their top and bottom edges. 

For a single plane, cut open along the upper and lower edges of the intervals 
(ui,vi) and (1*2,^2) there is a conformal mapping (whose explicit expression is not 
needed) to a cylinder of circumference W and length L. We introduce the modular 
parameter q = e - 2lTL / w which is a known function of the cross-ratio x of the original 
four points in the plane. The limit x — ¥ corresponds to q — ¥ 0. The upper and lower 
edges of the cuts are mapped into semicircular arcs (Six, Sir) and (S2L, S2_r) at each 
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Figure Al. The cut plane represented on the left can be mapped through a 
conformal transformation to a cylinder of length L and circumference W (right) 
and so modular parameter q = e~ 2wL / w , The upper and lower edges of the cuts 
are mapped into semicircular arcs (Six, S\r) and (S%l, S211) at each end of the 
open cylinder. 



end of the open cylinder. This is illustrated in Fig. |A1| We do this for each plane, 

labelled by j £ [l,n]. These are then connected cyclically by their edges so that 

is connected to S^ 1 ^ and S^l is connected to S^l • The case n — 4 is shown in 

Fig. go] 

Note, however, that this in general introduces new conical singularities at the 
points where the semicircles meet. The total angle subtended by a curve which 
encloses a singularity is now im, rather than 2im as in the original geometry. This 
changes the value of the trace anomaly at each singularity from (c/12)(n — (1/n)) to 
(c/12) ((n/2) — (2/n)J. The partition function in the coupled cylinder geometry thus 
has the form 

^cylinder ^ w ~(c/3) ((n/2)-(2/n)) j^yl ^ (A.l) 

However, the non-trivial dependence on q should be conformally invariant, that is 

F^(q(x))=F n (x), (A.2) 

where J- n {x) is the same function as in We assume they are both normalized so 
that 7£ yl (0) = ,F n (0) = 1. 

Note in particular that for n — 2 there is no trace anomaly. This because in this 
case, one of the cylinders can be reflected Lf>i?so that is now sewn onto S[ 2 ^ 
and is sewn onto S*^. For n = 2 we can also reflect the second cylinder in the 
plane of one of its ends. This gives a cylinder of length 2L with its opposite ends 
identified, that is a torus, with modular parameter e ~ i7rL / w — q 2 . As is well-known 
|32j the torus partition function encodes all the scaling dimensions of the theory, that 
is 

^Hl + E^T^- (A.3) 

fc^O 

This is consistent with our result above that, for n — 2, only terms involving the 
2-point functions of primary operators can arise in the small I expansion of J^, which 
has the form 

-F20C M^2/r 2 f A ^ + ---, (A.4) 

primaries 
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where dk — ^t-jj' an d the neglected terms are the contributions of descendant 

fields. Presumably the full identity of the above two expressions involves a complicated 
connection between these neglected terms and the higher order terms of the expansion 
of q(x) in powers of x. 

It is interesting to consider computing ^cylinder f or general n as a similar expansion 
in powers of q. This can done by considering the infinitesimal generator H of 
translations along each cylinder. As is well known [32], the eigenstates \k) of H 
are in 1-1 correspondence with the scaling fields of the CFT, with corresponding 
eigenvalue (27r /W)(Ak + Afe) — (irc/6W). The last term drops out if we normalize 
so that ^ yl (<? = 0) = 1. Now divide the Hilbert space into two subspaces Hl and 
T-Lr corresponding to the L and R halves of the cylinder. Then each state admits a 
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Schmidt decomposition 

|fc) = y^ y c ka \k;a) L ® \k;a) R . (A. 5) 

a 

Note the parity operator P which reflects L o R (and corresponds to complex 
conjugation in the planar geometry) gives an isomorphism between Hl and 7i R . 
However, the states \k;a)L and \k;a) R generated by the Schmidt procedure are not 
in general mapped into each other by this isomorphism. 
Since the states \k) form an orthonormal basis, 

^2ck a Ck' a ' L(k;a\k';a') L R (k; a\k'; a') R = 5 k k> • (A.6) 

For k = k' this follows from the orthonormality of the Schmidt states and the fact 
that *Yua c ka = 1> but m general it is not true that \k,a) and \k' ,a') are orthogonal 
for k 7^ k' . 

We make this decomposition for each end of each cylinder, labelled by j. The 
sewing procedure is then equivalent to taking inner products of states in the L space 
on cylinder j with those in the R space of cylinder j + 1. Explicitly 

^ yl (9)«E ( X II ' /■•..:<>,/••..::<>,., I g£* (A ** +S **\ (A.7) 

where we have dropped the distinction between 71 l and 7i R , assuming that these are 
identified using the isomorphism under P. However, when n is even we can relabel the 
basis for each even j so that \kj \ ctj) \kj\ ctj) and <H- (equivalent to reflecting 
cylinders with j even). The expression inside the modulus sign then becomes 

H II ^^{kj^jlkj+naj+i}] Yi c kj<x J {kj;a j \kj+i;a j+1 )\ . (A.8) 

{aj} \j odd J \j even j 

(For odd n we can do this with (n — l)/2 of the cylinders, leaving a single factor of 
(kj;ttj\kj+i;a j+ i).) 

For n = 2 we then see that the sum between the modulus sign has the same form 



as that in (A.6) and therefore we recover the torus partition function (A.3) 



It is illuminating to write (A.7) in another way. Introduce the operators 
Ml: >l,. 'Hl: 

M k = J2cka\k]a)(k;a\. (A.9) 



Then the expression in parentheses in (A.7) is simply 



3=1 

Note that under the parity operation M O M* . The fact that we can reflect an 
arbitrary number of the cylinders without changing the partition function means that 
we can, in the above matrix product, change M k . — > in any number of the factors, 
up to a sign. In fact the sign is determined by the parities of the states \kj): since 
P commutes with H, the eigenstates \k) can be chosen also to be eigenstates of P 
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with eigenvalues (— l) Afc j ~ Afc j . Note however that we should make the same set of 
replacements M — > M* for all {kj}. 

However, we can make a stronger statement, based on the observation that we can 
consider the modulus qj of the jth cylinder to depend on j, and also to be complex, 
corresponding to twisting the cylinder before the sewing procedure. This no longer 
corresponds to computing the Reny i entropy, but still makes mathematical sense. In 

is generalized to the form JT qf'qf 3 . Therefore 



this case the g-dependence in (A. 7 



the identity between the different expressions obtained by reflecting each cylinder 
independently must also hold for each set {kj} with the same values of {A*,. } and 
{Afc.}. The consequence of this is most easily seen in the case when all the {kj} 
correspond to primary operators. Then, barring accidental degeneracies, 

Tr (M fel M fc3 ---M fen ) , (A.ll) 

is invariant, up to a sign determined by the overall parity, on replacing any subset 
of the M kj by their hermitian conjugates. Since this must be true for arbitrary n, 
the simplest solution is that M\ — ±M k . In the more general case, there is a finite- 
dimensional space of descendent operators with the same values of (Afc, A&), and we 
have 



Ml=YA kk ,M k ,, (A. 12) 



k 



where the sum is restricted to operators within this space, and A* A = 1. 

We can express the results of Sec.[3]in this formalism. The case n = 1 corresponds 
to a sphere, which has no modulus and therefore there should be no (/-dependence. 
This implies that TrM k cx 5 k o- Similarly, for n — 2, we have Tr M k M k i cx 5 k k'- 

If we take all the kj — we get the leading term in the partition function as 
q^O: 

TrM™ oc W^ (c/3) ( ( " /2) - (2/n) ) . (A.13) 
For even n, we can also write the left side as 

Tt{M Q MlY' 2 = Y,<%a- (A- 14) 

OL 

For odd n we get Tr (Cq^ 1 M ) , where Co is a diagonal matrix with entries co a . This 
gives information about the spectrum of Mo and, since this is valid for all n, suggests 
that its non-zero eigenvalues are simply cofc. However, since Hl is infinite-dimensional, 
this may not hold in a strong sense. 

Once this overall W dependence is factored out, we see, that, at least for primary 
operators 

n 

Tr J] M kj cx (J] fcj (e 2 ^>)) c , (A. 15) 

i=i j 

where the rhs is the n-point function at the nth roots of unity in the plane. It would 
be interesting to explore the properties of the operators M k further, as they appear 
to encode the entanglement properties of the theory in a slightly different way. The 



operators (p k . in Eq. (A. 15) can be thought as "propagating" in the cylinders as 
pictorially shown in Fig. A*2 (where they have been called Oj). 
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